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Abstract. Resonance production and decay into pion pairs is simulated in a non-extensive quark matter witii 
multi-particle interactions. Final state pion spectra are found to take the form of the Tsallis distribution, in ac- 
cordance with measurements. It has also been shown that, if a large number of particles with these multi-particle 
interactions are constrained to a constant energy hyper-surface in phase space, the one-particle distribution is the 
Tsallis distribution. 



1 Introduction 

Transverse hadron spectra measured in high energy colli- 
sions in the last three decades, fit to the Tsallis distribution 
(TS) (see Refs. im- lfT2l for proton-proton (pp), proton- 
antiproton (pp) and nucleus-nucleus (AA) collisions and 
Refs. 1 13 14| for e^e" collisions). On the theoretical side, 
there are many proposals on the emergence of the TS dis- 
tribution. In kinetic theory, the collision term of the Boltz- 
mann-equation ifTSl - lfTTll . or the noise term of the Langevin- 
equation | lY] can be generalised in a way, in which the 
TS distribution is the stationary soUution. In equilibrium 
thermodynamics, the Maximum Entropy Principle (MEP) 
together with a generalisation of the Shannon-entropy for- 
mula (the Tsallis-entropy, see Ref. 1 19 1) also lead to the TS 
distribution, as a generalisation of the canonical Boltzmann- 
Gibbs distribution (BG). The TS distribution can also be 
derived from the MEP by introducing special interactions, 
while, leaving the original Shannon-entropy unaltered II II20I . 

In all the above cases, equilibrium, or at least stationar- 
ity is assumed, while the question of equilibration in high 
energy collisions is still a subject of intense debate. How- 
ever, examining the mathematical foundations of statistical 
physics |21 22 1, it turns out, that statistical physical distri- 
butions may be used not only in equilibrium. Whenever 
a system is composed of idependent and identically dis- 
tributed particles, and the total energy of the system is con- 
served, the one-particle distribution may be approximated 
by the canonical distribution in the limit of a large number 
of particles. The reason for this is that Central Limit The- 
orems (CLT) do not deal with the issue, whether particles 
have the same distribution, because they are thermalised, 
or because they are produced via the same process. 

In this paper, we discuss resonance production in a 
non-extensive quark-gluon medium, in which multiparticle 



interactions of the form of Eq. Q are present. In Sect. [2] 
we show, that (apart from a phase-space factor) the one- 
particle energy distribution in such a medium is the TS 
distribution in the limit of a large number of particles. In 
the calculation we use only probability theory, thus the re- 
sult is valid not only in equilibrium. We also outline some 
advantages and disadvantages of the probability theory ap- 
proach to the thermodynamical one. 

In Sect. |3] we introduce a model in which a non-exten- 
sive quark matter (QM) evolves in time via the collisions of 
randomly chosen quark pairs. In the collisions, the 3-mo- 
mentum, and the total (non-extensive) energy of the sys- 
tem is conserved. In Sect. 3.1 we present a way, in which 
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hadron resonances may leave the QM and decay into pion 
pairs without the violation of the conservation of the total 
energy and momentum of the system. 

Sect. |4] contains the resulting final-state pion spectra 
and the mass spectrum of the resonances. Sect. |5] contains 
our concluding remarks. 



2 Probability Theory, Maximum Entropy 
Principle and the Tsallis Distribution 

In this section we review the derivation of the canonical 
BG distribution using probability theory. During the de- 
duction, we allow for a constraint on the N particle phase 
space that is more general than the conservation of the sum 
of the one-particle energies. This way, we may obtain a 
wider class of distributions covering those having power- 
law asymptotics as well. 

Let us assme, that a system is composed of a large 
number of particles, all being independent and identically 
distributed (iid.). For example let the particles have the 
same phase space structure with momentum-space distri- 
bution, t/F(pi) (in a homogenious, isotropic ensemble 
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dF{pi) oc d^pi). In this case the probabiHty that a parti- 
cle has energy, e, while the total energy of the system, E is 
fixed, is 



dF{p) Qn{E) 

with E' being the energy of - 1 particles (when particles 
do not interact, E' - E - e) and QNiE) being the phase 
space volume of particles restricted to the constant en- 
ergy hyper-surface: 



Note, that /3 and thus ;F also depend on L{E)/N. Applying 
Eq. (j7]i, the one-particle energy distribution (Eq. ([T]i) in the 
large N limit, is approximately 



dVN 
dF(p) 



exp { -{N -\)T 



HE') 



N-l 



L(E) 



N 



(10) 

with L{E') - L{E) - L{e). Exploiting that £ » e, in the 
first term in the bracket, we may Taylor-expand T around 
LiE)/N. It follows from Eq. dSli that 



via the constraint. 



/n 



dF(pi) 6 



Y^L(ej)-L(E) 



(2) 



(3) 



For example, if the particles are non-interacting, 2 ej - E, 
thus L(e) - e. In the model, reported in IfTTI . the following 
multiparticle interaction has been proposed: 



Em - e\ + ■ ■ ■ + €n 

+ a{e]_e2 + ... + eN-\eN) 



+ a ei ■ ■ ■ e^. 



(4) 



The above rather complicated energy formula can be turned 
into a simlpe addition of the form of Eq. (|3]l by a so-called 
formal logarithm Il20l . 



(5) 



Lie) = -ln(l +ae). 
a 



Moreover, the constrained particle phase space volume, 
Eq. (|2| can be factorised using In 6(x) - J ds exp(i s x): 



oo 

QAE) ^ ^ ds exp|-A^ 



• ^^^^ u ^ 



(6) 



with the logarithmic generator function of the one-particle 
phasespace distribution, t(s) - In ^ J (iF(p) exp [-/ s L(e)]). 
In the large limit, Eq. (j6]) may be approximated by the 
saddle-point method: 



QN(E)^tx^[-Nr[s*]] 
where T is the free energy per particle, 

r{s*]^-is*^-i{s*), 

and the inverse temperature /3 = i s* minimises 'F: 

LjE) _ JdF(p)L(e)exp[-(5L(€)] 
N " JdF(p)exp[-pL(e)] 



(7) 



(8) 



(9) 



r' 



LjE) 
N 



^-is- +\-i^-i'is*) \ s*'^-fi. (11) 



The expression in the bracket in Eq. ( [TT| vanishes because 
of Eq. (|9|. Consequently, if we neglect the terms in the 
Taylor-series that are proportional to 1 jN (or smaller), the 



one-particle distribution, Eq. ( 10 1 gives the distribution 



dpK 



exp{-^L(e)) 



dF{p) J dF(p) exp{-/3L(e)} 



(12) 



The above calculation is similar to the variational met- 
hod based on the Maximum Entropy Principle. There, a 
free energy functional 'P[f] is constructed from the one- 
particle distribution /, that contains an entropy term, S [/] 
and constraints C[f]: 



/3 <!>[/] = -S[f]+/3C[f]. 



(13) 



There are several proposed formulas on S [/], however there 
is no recipe on how to construct the appropriate S [/] for 
a system defined by a given Hamiltonian. The following 
choices of C[f] and S [/] both lead to the TS distribution: 



S[f] 



j dF(p)f{e)ln[f(e)] 
C[f] = J dF(p)L(e)f(€)- ^^^^ 



N 



(14) 



S[f] = 



1- JdF(p)f%6) 



C[f] 



-f 



q- 



1 



dF{p)ef(e)- 



N 



(15) 



In Eq. ([T4]i, the entropy is additive for independent par- 
ticles, while interactions of the form of Eq. Q are pres- 
sumed, thus the mean one-particle energy is constrained 
through the formal logarithm. 

In Eq. (15 I, particles are assumed to be non-interacting, 
thus the mean energy per particle is set to be E /N, however, 
a generalised entropy formula accounts for correlations. 

Both the Maximum Entropy Principle and the probabil- 
ity theoretical approach are based on finding the minimum 
of a free energy, however, there are two main differences 
between them: 
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Fig. 1. The E] - E2 plane, in case, when the total momentum of 
the incoming quarks P is sufficiantly small. Because of the trian- 
gle inequality among pi, p2 and P, the energies of the outgoing 
quarks (Ei and £2) have to be chosen from the filled area. Fur- 
thermore, E[ and E2 have to satisfy Eq. ^19\ (solid red curve). 

1) Behind the Maximum Entropy Principle, there is the 
assumption of thermal equilibrium, while in the prob- 
ability theoretical approach, there is no need for such 
an assumption; 

2) In the probability theoretical approach, particles need 
to be independent, while this is not required in a model 
based on the Maximum Entropy Principle. 



3 Parton Collision Cascade Simulation and 
Pion Production 

In this section, we use a model reported in f\T\. We con- 
sider massless quarks (e,- = |pi|) with interactions of the 
form of Eq. Q. The ensemble has an initial momentum 
distribution that is homogenious inside a Fermi sphere of 
radius pp. The momentum distribution evolves in time via 
collisions of randomly chosen pairs. In pair collisions the 
3-momentum and the formal logarithm of the energy are 
conserved: 



in 



Fig. 2. The £1 - £2 plane, in case, when the total momentum of 
the incoming quarks P is sufficiantly large. Because of the trian- 
gle inequality among pi, p2 and P, the energies of the outgoing 
quarks (£1 and £2) have to be chosen from the filled area. Fur- 
thermore, El and £2 have to satisfy Eq. \\9\ (solid red curve). 



with the total momenta and energy of the incoming quarks: 
P and Ei„. Integrating out for p2, we obtain the energy dis- 
tribution of the first outgoing quark: 



p{E,) 



E\ (Ej,, - El) 

(1 +fl£l)2 



(18) 



Once El is obtained, the energy of the second outgoing 
quark is 



El = 



E ill E 



1 +aE 



(19) 



Pi + P2 = P3 + P4 , 

L(ei) + L(e2) = L(e3) + L(e4). 



(16) 



(The second row in Eq. ( 16 1 also ensures the conservation 



of not just the pair energies but that of the total system.) 
Thus the momenta of the outgoing quarks are chosen ran- 
domly according to the distribution 



The paiticles are massless, thus the equation pi -H p2 = P 
defines the angles between their momenta and P. However, 
£1 may not take any value between and Ei„, because pi, 
P2 and P must obey the triangle inequality. 

• If P < I ^ VI + aEi„ - 1) , we have to choose £1 from 

the interval £, e ±f - i + V(f f + ^ + V (^^^^^ 

• If P > VI + Ei„ - 1) , £■] must satisfy the con- 



ditions: £1 e 

Eii^± 



(Fig. 



21, 



+ -k + — , and that 



dw - d^pi d^p2 S (pi -H P2 - P) 

X 5{Ei + E2 + aEiEi - £■,„) 



(17) 



An interesting feature of this model is its exponential 
non-extensivity. Let us estimate the dependence of the total 
energy of the system on the number of its constituents in 
the case, when the mean energy per particle, e is fixed. 
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Fig. 3. Schematic picture of resonance production in a non- 
extensive medium. Spots represent quarks, curves represent in- 
teractions among them. 

From Eq. ([3]l, the total enrgy of the system estimated from 
the equipartition principle, e, - e, is 



E ^ L-'[NL(e)] 



1 



-(1 +aey 



(20) 



3.1 Resonance Production in a Non-extensive 
Medium 
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Eq. (20 1 illustrates that, when a color-neutral quark anti- 



quark pair leaves the system, and forms a hadron reso- 
nance, it carries a big amount of interaction energy, (see 
Fig.[3]as a schematic picture of resonance production) Be- 
cause of the 3 -momentum conservation, this hadron can 
not be on-shell. For simplicity, we do not distinguish be- 
tween quark flavours, and consider only a pion pair pro- 
duction in the decay of the resonance. Since the resonance 
{R) has left the quark matter (QM), R does not interact with 
the QM any more. Thus the energy of R is simply the dif- 
ference of the total energy of the QM before (En) and after 
{Eff_2) the formation of R: 



Eb — En — Em-?. — El 



1 +aE 



N 



' 1 + aEi2 



(21) 



with £12 - Ei+E2+aE]E2 being the energy of the qq pair 
without its interaction with the rest of the QM. The second 
equality and the last term in Eq. pT) follow from Eq. Q. 
This shows, that the energy of R can be much greater than 
that of the qq pair, while the momentum of R is exactly the 
same of that of the qq pair Consequently, the production 
of big resonance masses is allowed, as shown in Fig. |4] 
The rest mass of R then contributes to the kinetic energies 
of the pions when R decays, resulting in a long tail of the 
momentum distribution of the pions, as shown in Figs.|5]|7] 



4 Results 

The model described in Sect. [5] converges to the TS distri- 
bution (Eq. ( 12 1 in the case, when L(e) - {I /a) ln(l -1- ae)) 



dp 



- Ap (l + a£) 



(22) 



when started from a homogeneously filled Fermi sphere. 



Fig. 4. Simulated resonance mass distribution from the model de- 
scribed in Sect. [3] The histogram in the figure is the sum of statis- 
tics collected in 1000 events with 50 resonance produced in each 
event. Through out the simulation, a = I was used. 



to the model as well, the results shown in Figs. |4||7] are ob- 
tained. As a result of the large interaction energy caused 
by the multi-particle interaction therms in Eq. Q, the pro- 
duced resonances may have large masses, as can be seen in 
Fig.|4] The mass distribution of the resonances, dNnldMn 
shows power-law behaviour for two orders of magnitude 
in the mass Mr. Low masses, Mr < 2m„ are not allowed 
throughout the formation. 

When the resonances decay, their rest energy contributes 
the kinetic enrgies of the produced pions. This causes long 
power-law tails of the final state pion spectra shown in 
Figs. 5]|7 Fig. [5] shows the dependence of the final state 
pion spectrum on the number of quarks (which is equal to 
the number of final state pions in the process qq ^ R 
TT^n^) in the QM, when the interaction measure, a - I. 



Apparently, the spectra take the form of Eq. ( 22 1 in the 



canonical region, where the one-particle energy is much 
smaller, than the total energy of the system, e Eff. For 
larger energies, e ^ eo, however the spectra have a cut. eo 
is proportional to the total energy, and grows like eo ~ e'^. 



Figs. |6] and |7] show the dependence of the spectrum of 
events with fixed multiplicity (A^ - 2000) on the interac- 
tion measure a. Apparently, the spectra change from BG 
to TS, as the interaction measure grows from a - 0.01 to 
a = 0.2. The spectra with a = 0.01 and a - 0.02 are nearly 
BG distributions (apart from the cuts above eo), while spec- 
tra with a = 0.05, a = 0.1 and a - 0.2 are well-developed 



When resonance production (described in Sect. 3. 1 1 is added TS distributions (Eq. ( 22 1) 
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Fig. 5. Simulated + n" spactra from the model described in 
Sect. |3] The colour encoding distinguishes between the spectra 
of events with different multiplicities. Each histogram contains 
statistics of 1000 events. Through out the simulations, a = 1 was 
used. 



Fig. 7. Simulated + n" spactra from the model described in 
Sect. [3] shown in a log-log plot. The colour encoding distin- 
guishes between the spectra of events with different interaction 
measures, a. Each histogram contains statistics of 100 events with 
2000 ;r-s per event. 
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Fig. 6. Simulated + tC spactra from the model described in 
Sect.|3] shown in a log-lin plot. The colour encoding distinguishes 
between the spectra of events with different interaction measures, 
a. Each histogram contains statistics of 100 events with 2000 ff-s 
per event. 



5 Conclusions 

In this paper, we discussed hadronisation in a non-extensive 
quark matter with multipaiticle interactions of the form of 
Eq. (|4]). In Sect.[2]we have shown that the one-particle dis- 
tribution is the Tsallis distribution in a system, which is 
composed of independent and identically distributed parti- 
cles, that are constrained on a constant energy hyper-surface 
in phase space by Eq. (j4]). (We note, that this statement is 
valid only in the limit of a large number of particles, and 
if the one-particle energy is much smaller than the total 
energy of the system.) Furthermore, this statement holds 
for non-equilibrium systems as well (this, being its advan- 
tage over the approach based on the Maximum Entropy 
Principle (MEP)), however, the MEP may be used for non- 
independent particles too. 

In Sect. [3] we have presented a hadronisation process, 
that respects the conservation of the total 3-momentum and 
non-extensive energy of the system. In the process, first 
resonances are formed, then decay into pion pairs. We have 
found that the resonances may have large masses, due to 
the large interaction energies in the quark matter, caused 
by the multiparticle interactions. This leads to the power- 
law tailed resonance mass distribution shown in Fig.|4] The 
resonance masses then contribute to the kinetic energies of 
the pions, resulting in long, power-law tailed pion spectra 
presented in Sect. |4] 

Since both our simulated and the measures pion spectra 
take the form of the Tsallis distribution, we may conclude, 
that the parameters of our model (mean energy per parti- 
cle and interaction measure (a)) may be tuned so that our 
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results and experimental data be in accordance with each 
other. 
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